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1.0  SUMMARY 


The  failure  of  Reaction  Wheels  (RWs)  in  an  array  can  impair  the  spacecraft's  ability  to  perform 
imaging  missions,  during  which  a  prescribed  inertial  orientation  has  to  be  achieved  and  maintained 
with  high  accuracy  and  precision.  Even  though  spacecraft  are  equipped  with  four  or  more  RWs  for 
redundancy,  the  possibility  of  multiple  RW  failures  does  exist.  Even  though  the  spacecraft  may 
still  retain  thruster-based  attitude  control  capability,  the  attitude  control  accuracy  that  thrusters  are 
able  to  attain  and  maintain  is  typically  insufficient.  In  addition,  the  use  of  thrusters  involves  using 
fuel,  which  shortens  the  spacecraft's  life.  Hence  with  RW  failures,  the  spacecraft  becomes 
underactuated.  Recent  missions,  such  as  Kepler  and  FEISE  (Ear  Ultraviolet  Spectroscopic 
Explorer),  are  exemplar  of  these  situations  [1,2].  The  focus  of  this  research  was  on  the 
development  of  several  approaches  that  could  be  applied  to  an  underactuated  spacecraft,  satisfying 
high  precision  pointing  requirements  and  compensating  for  disturbances.  Three  methods  were 
explored; 


1.  Utilizing  Solar  Radiation  Pressure  (SRP)  to  restore  linear  controllability  to  a 
spacecraft  with  two  RWs. 

2.  Creating  a  set-optimization  technique  that  computes  drift  sets.  These  sets  represent 
all  initial  orientations  such  that,  if  the  spacecraft  starts  in  the  set,  it  would  satisfy 
constraints  up  to  a  given  time  regardless  of  drift. 

3.  Developing  hybrid  and  Model  Predictive  Controller  (MPC)  approaches  for 
reorientation  and  disturbance  compensation  in  underactuated  spacecraft. 

Multiple  publications  have  been  published,  written,  or  are  in  the  process  of  being  written,  and  are 
related  to  the  subject  matter  of  this  report.  All  such  publications  are  listed  in  the  Bibliography. 


2.0  INTRODUCTION 


In  the  case  of  an  underactuated  spacecraft  with  two  internal  torque  actuators,  the  linearized 
dynamics  are  not  controllable.  In  [3]  the  nonlinear  system  is  said  not  to  be  accessible  with  two  or 
fewer  momentum  wheels  due  to  the  conservation  of  angular  momentum.  In  many  studies,  the 
angular  momentum  of  the  spacecraft  system  is  assumed  to  be  zero.  In  [4],  under  the  zero  angular 
momentum  assumption,  two  nonlinear  control  techniques  are  presented.  The  first  is  a  finite-time, 
discontinuous  feedback  law  that  induces  a  sequence  of  rotations.  The  second  involves  a 
diffeomorphic  transformation  to  a  simpler  set  of  equations  of  motion  that  facilitates  controller 
design.  These  control  strategies  are  discontinuous,  as  the  attitude  of  an  underactuated  spacecraft 
cannot  be  smoothly  or  continuously  stabilized  by  any  time-invariant  feedback  law  [5,6,7].  Control 
techniques  using  the  zero  angular  momentum  assumption,  however,  are  of  limited  practical  use. 
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Under  this  assumption,  a  spacecraft  with  less  than  four  RWs  must  spin  down  the  RW  speeds  to 
zero,  which  as  a  consequence  decreases  accuracy  and  may  further  decrease  their  operational  life. 


3.0  METHODS,  ASSUMPTIONS,  AND  PROCEDURES 

3.1  Assumptions 

For  this  work,  a  spacecraft  configuration  consisting  of  a  spacecraft  bus  and  four  RWs  is 
considered.  During  the  spacecraft's  mission,  two  of  the  RWs  will  fail,  leaving  two  operational 
RWs  about  the  first  two  principal  axes  of  the  spacecraft  bus.  The  RWs  are  assumed  to  be  thin  and 
symmetric.  The  orientation  of  the  spacecraft  is  characterized  by  3-2-1  Euler  Angles  yaw  (\|/),  pitch 
(0),  and  roll  (cp),  with  the  kinematics  prescribed  by 


0  =  M(0)a),  (1) 

where0  is  a  vector  of  Euler  angles,  o)  is  a  vector  of  angular  velocities,  and  M(0)  is  a  matrix  which 
is  a  function  of0.  Any  three  parameter  parameterization  of  orientation  has  singularities.  It  is 
assumed  that  the  spacecraft's  mission  does  not  require  large  attitude  transients  and  never 
approaches  singularity. 

The  dynamics  of  the  spacecraft  system  are  governed  by  Euler's  equation, 

(J  +UWW^)(b  =  -ajx{(J  +J^WW'^)aj  +  UWv)  -  J^Wv  +  (2) 

where  /is  inertia  matrix  of  the  spacecraft  bus,  /^^is  the  inertia  of  the  RWs  about  the  spin  axis,  W  is 
the  orientation  matrix  of  the  RW  spin  axes  with  respect  to  the  spacecraft  bus,v  is  the  vector  of  the 
reaction  wheel  speeds,  and  Mg^f-is  the  vector  of  external  moments  acting  on  the  system.  Eor  all  the 
methods  discussed,  the  external  moments  are  limited  to  control  moments  and  SRP  torques.  In  the 
development  of  the  MPC  approach,  it  will  be  assumed  that  there  are  no  external  moments  acting 
on  the  spacecraft. 

Specifically  for  the  MPC  approach,  the  8-dimensional  state  vector  for  a  spacecraft  with  two  RWs 
can  be  reduced  to  a  5 -dimensional  set  of  equations  due  to  the  law  of  angular  momentum 
conservation.  These  equations,  called  the  reduced  attitude  equations,  are 


0  =  M(0)a  +J^WW^rKOBGH  -J^Wv),  (3) 

whereCgc  is  the  direction  cosine  matrix  from  the  inertial  frame  to  the  spacecraft  frame  and  H  is  the 
inertial  angular  momentum  vector. 

The  SRP  model  used  in  this  work  is  based  upon  that  of  [8].  Since  the  distance  from  the  spacecraft 
to  the  Sun  is  large,  SRP  is  assumed  to  act  evenly  across  the  same  panel.  The  total  SRP  exerted  on 
a  panel;  is 
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Pj  =  -a{Unj  ■  Us){Unj  +  PjUs), 


(4) 


wherea  and  Pj  are  constants  which  are  functions  of  the  distance  from  the  Sun,  solar  flux,  speed  of 
light,  and  diffusion  coefficient  of  the  panel,  u^j  is  the  normal  to  panel  j  pointed  outward,  and  Ug  is 
the  unit  vector  from  the  center-of-mass  to  the  Sun.  The  solar  radiation  torque  exerted  by  panel  j  is 
given  by 


^jsrp  ~  iPjO  ^co)  ^  (^) 

whereryo  is  the  vector  from  the  center-of-mass  to  the  center  of  panel  j,  Tqq  is  the  vector  from  the 
center-of-pressure  to  the  center-of-mass,  and  Aj  is  the  area  of  panel  j.  Since  SRP  is  additive  across 
all  faces,  the  total  SRP  torque  exerted  on  the  spacecraft  is  the  sum  of  the  torques  due  to  all 
individual  panels. 


3,2  Methods 

3,2.1  Recovering  Linear  Controllability  with  Solar  Radiation  Pressure 

For  this  approach.  Solar  Radiation  Pressure  (SRP)  torques  are  included  in  the  dynamics  of  the 
spacecraft  model.  For  small  deviation  from  inertial  pointing,  the  linearized  model  can  be  used 
effectively  for  analysis  and  controller  design.  While  the  open-loop  system  is  unstable,  our  results 
show  that  under  appropriate  assumptions,  which  are  outlined  in  [9,10],  the  linear  system  becomes 
controllable.  The  assumptions  include  the  solar  radiation  pressure  center  being  different  from  the 
center-of-mass  of  the  spacecraft. 

With  linear  controllability  regained,  spacecraft  stabilization  can  be  achieved  by  conventional 
control  schemes.  A  Linear  Quadratic  (LQ)  approach  was  first  applied  due  to  its  robustness,  its 
optimal  control  properties,  and  its  familiarity  to  aerospace  engineers.  A  pole  placement  scheme 
was  also  used  to  improve  convergence  time,  since  with  LQ  the  closed-loop  response  was  relatively 
slow.  Overall,  by  taking  advantage  of  the  change  in  the  dynamics  induced  by  SRP  torques, 
two  RWs  are  able  to  slowly  correct  the  attitude  errors  over  time. 


After  [10]  was  published,  the  press  release  [11]  suggested  that  SRP  is  actually  used  (in  an 
unspecified  control  scheme)  to  restore  Kepler's  mission  controllability.  The  controllability 
analysis  and  results  from  this  work,  obtained  independently  of  [11],  are  thus  indirectly 
corroborated  by  experimental  evidence  in  [1 1]. 
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3.2,2  Disturbance  Torque  Counteraction  with  RWs  and  SRP  Torques 


In  this  part  of  the  research  we  focused  on  characterizing  the  ability  of  the  spacecraft  to  perform 
inertial  pointing  while  compensating  for  disturbance  torques  caused  by  SRP.  This  is  done  by  the 
computation  of  drift  sets,  which  characterize  the  set  of  all  initial  conditions  such  that  if  the 
spacecraft  starts  in  the  set,  it  will  achieve  the  desired  inertial  pointing  given  constraints  on  pointing 
accuracy,  pointing  time,  angular  velocity,  RW  speeds,  and  RW  accelerations.  Since  small 
deviations  are  considered,  a  linear  discrete -time  model  is  derived,  with  the  form 

x(k  +  1)  =  Ax{k)  +  Bu{k')  +  w(/c),  (6) 

where  x  is  the  state  (vector  of  attitude  Euler  angles  errors  and  angular  velocity  errors),  u  E  U  is  the 
control  (vector  of  two  RW's  torques)  and  w  EW  represents  the  effects  of  unmeasured 
disturbances/uncertainties  bounded  by  the  set  W .  The  constraints  on  the  system  are  represented  by 
n  linear  equality  constraints  on  the  state, 

HiX{k)  <  hit  =  1, ...  ,n.  (7) 

These  linear  inequality  constraints  can  be  used  to  construct  the  setO/^,  which  is  defined  as 

Ok  =  {^0  S  R^\HiX{k)  <  hi,i  =  1, ...  ,n} ,  (8) 

where  Xq  —  x(0).  The  set  0^.represents  the  set  of  all  initial  condition  such  that  the  spacecraft  will 
satisfy  constraints  at  time  k.  The  intersection  of  all  sets  Ok  up  to  a  given  discrete  time  instant  p 
will  give  the  set  of  all  initial  conditions  such  that  the  spacecraft  will  not  violate  constraints  up  to 
and  including  the  time  instant,  p.  In  addition,  these  sets  can  be  easily  represented  and  visualized 
by  polytopes.  Future  work  will  involve  using  these  drift  sets  in  order  to  construct  a  drift  reaction 
controller. 


3,2.3  Hybrid  and  Model  Predictive  Controller  for  Reorientation  and  Disturbance 
Compensation  in  Underactuated  Spacecraft 

The  hybrid  approach  developed  for  the  underactuated  spacecraft  problem  is  based  on  the  feedback 
stabilization  techniques  developed  in  [12]  and  [13].  These  techniques  exploit  the  decomposition  of 
the  system  variables  into  base  variables  and  fiber  variables.  The  base  variables  are  stabilized  to 
periodic  motions  with  feedback,  and  the  parameters  of  these  periodic  motions  are  adjusted  at 
discrete  time  instants  to  induce  the  evolution  of  the  fiber  variables  towards  the  desired  equilibrium. 
In  the  case  when  the  spacecraft  is  actuated  by  two  external  moments,  the  attitude  and  angular 
velocity  of  the  two  actuated  axes  are  treated  as  base  variables  and  the  attitude  and  angular  velocity 
in  the  third,  uncontrolled  axis  are  treated  as  fiber  variables.  In  the  case  when  the  spacecraft  is 
actuated  by  two  internal  moments  provided  by  RWs,  the  base  variables  are  the  same  as  in  the  two 
external  moment  case,  but  the  only  fiber  variable  considered  is  the  attitude  about  the  uncontrolled 
axis.  When  the  system  was  controlled  by  two  external  moments,  the  system  could  reject  SRP. 
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When  internal  moment  control  was  used,  a  periodic  reset  scheme  was  implemented  for  a  controlled 
oscillation  about  the  uncontrollable  axis. 

For  the  MFC  approach,  the  reduced  attitude  equations  were  used,  and  hence  it  was  assumed  there 
were  no  external  moments  acting  on  the  system.  It  was  proven  that  for  any  nonzero  angular 
momentum,  the  reduced  attitude  equations  are  linearly  controllable.  An  MFC  control  law  can  then 
be  applied  to  stabilize  the  desired  equilibrium  as  long  as  the  angular  momentum  vector  expressed 
in  the  spacecraft  bus  fixed  frame  is  in  the  range  of  the  matrix  W.  If  the  angular  momentum  is  not 
in  the  range,  then  under  suitable  conditions,  the  spacecraft  will  converge  to  a  nonzero  equilibrium. 

Both  the  hybrid  and  MFC  method  above  work  for  nonzero  angular  momentum.  The  topic  of 
underactuated  spacecraft  control  when  the  angular  momentum  is  non-zero  has  not  been  studied 
extensively  in  the  literature,  but  is  of  great  importance.  By  creating  controllers  that  include  and 
utilize  nonzero  angular  momentum,  the  operational  life  of  the  RWs  can  be  extended.  Studying 
nonzero  angular  momentum  also  increases  the  scope  of  the  spacecraft  underactuated  problem,  as 
zero  angular  momentum  is  a  limiting  case. 

For  zero  angular  momentum,  the  system  is  linearly  uncontrollable.  In  fact,  the  attitude  problem 
with  two  RWs  and  zero  angular  momentum  cannot  be  stabilized  through  any  smooth  time- 
invariant  feedback  law.  This  is  a  consequence  of  Brocketfs  theorem  in  continuous  time  [5],  which 
can  be  extended  to  discrete  time  [14].  However,  nonlinear  MFC  can  result  in  a  discontinuous 
control  trajectory.  Therefore,  a  nonlinear  solver  will  be  applied  in  future  work. 


4.0  RESULTS  and  DISCUSSIONS 

The  simulations  presented  use  the  spacecraft  parameters  listed  in  Table  1 .  In  Table  1,  O^^^jis  the 
solar  flux,  Qt//  is  the  diffusion  coefficient,  assumed  to  be  the  same  for  all  panels,  and  c  is  the 
speed  of  light.  It  is  assumed  that  the  two  RWs  remaining  after  failure  are  those  two  whose  spin 
axis  orientation  correspond  to  the  first  two  columns  ofW,  (i.e.  those  aligned  with  two  of  the  three 
principal  moments  of  inertia). 


Table  1.  Spacecraft  Parameters 


Parameter 

Units 

Value 

J 

kg  iiF 

tliag(430,  1210.  1300) 

Jw 

kg  nr^ 

0.043 

4. 

m2 

12.5 

Ay 

nF 

10 

A, 

nF 

5 

'10  0  l/\/3' 

W 

- 

0  1  0  1/C3 

0  0  1  l/\/3_ 

*^sun 

W/nF 

1367 

Cdiff 

- 

0.2 

C 

m/sec 

299792458.0 
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4,1  Recovering  Linear  Controllability  with  Solar  Radiation  Pressure 

The  following  simulation  illustrate  the  case  when  RW  3  fails  first,  followed  by  RW  4.  The  center- 
of-pressure  is  offset  from  the  center-of-mass  by  0.5  meters.  The  responses  are  shown  in  Figure  1. 
There  is  a  much  larger  disturbance  to  orientation  when  the  second  wheel  fails.  The  controller 
manages  the  first  wheel  failure  quickly,  then  reconfigures  and  handles  the  second  wheel  failure 
over  a  longer  period  of  time  coordinating  two  operational  RWs  in  presence  of  SRP. 
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(a)  (b) 

Figure  1,  (a)  Euler  Angles  and  (b)  RW  Speeds  for  When  2  RWs  Fail  and  SRP  is  Used  to 

Stabilize  the  Spacecraft 

4,2  Disturbance  Torque  Counteraction  with  RWs  and  SRP  Torques 


In  many  spacecraft  missions,  it  is  desirable  to  perform  rest-to-rest  maneuvers.  By  setting  the 
angular  velocity  entries  of  Xq  to  zero,  the  linear  equalities  become  functions  of  only  orientation. 
Therefore  the  drift  sets  can  be  easily  visualized.  See  Figures  2  and  3.  For  these  drift  sets,  the 
limits  on  spacecraft  orientation  and  angular  velocity  are  1  arc  second  and  0.0052  rad/sec 
respectively.  Since  an  linear  quadratic  regulator  controller  is  used,  the  RW  accelerations  can  be 
directly  obtained  from  the  orientation  and  angular  velocity.  Therefore,  a  limit  of  10  rad/sec  on 
RW  acceleration  is  imposed.  In  both  figures,  the  red  box  represents  the  constraints  on  pointing 
while  the  blue  box  is  the  set  of  all  initial  conditions  the  spacecraft  can  start  and  not  violate 
constraints  after  k  time  steps.  As  k  increases,  the  blue  box  will  decrease  and  either  converge  to  a 
set,  as  in  Figure  2,  or  vanish  if  constraints  are  too  stringent,  as  in  Figure  3. 
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Figure  2,  Drift  Sets  Calculation  Over  (a)  10  min  (b)  20  min  (c)  30  min  (d)  40  min  When  the 
Controller  Can  Fully  Compensate  for  SRP  Torques 


Figure  3,  Drift  Sets  Calculation  Over  (a)  2  min  (b)  4  min  (c)  6  min  (d)  8  min  When  the 
Controller  Cannot  Fully  Compensate  for  SRP  Torques 
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4,3,  Hybrid  and  Model  Predictive  Controller  for  Reorientation  and  Disturbance 
Compensation  in  Underactuated  Spacecraft 

Figure  4  illustrates  the  use  of  the  hybrid  seheme  to  stabilize  to  the  origin  when  there  is  an  offset  in 
the  uneontrollable  Euler  angle.  The  initial  eonditions  are  sueh  that  the  angular  momentum  is 
nonzero.  As  the  fiber  variable  deereases  to  zero,  the  amplitude  of  the  base  variables  deerease  and 
they  go  to  zero  as  well,  eonverging  to  the  eorrect  orientation. 


(a)  (b) 

Figure  4,  Response  of  Spacecraft  Controlled  by  two  Internal  Moments  Produced  by  RWs 
and  Produced  by  tbe  Hybrid  Controller  to  (a)  Euler  Angles,  (e)  Excitation  Magnitude 


Figure  5  demonstrates  the  MPC  solution  when  applied  to  the  redueed  attitude  equations.  The 
advantage  to  MPC  is  that  it  allows  for  eonstraint  enforeement.  Therefore,  constraints  are  added  to 
prevent  zero  crossings  of  RW  speeds  and  limits  of  RW  accelerations.  In  this  simulation,  there  is  5 
rad/sec^limit  of  RW  acceleration.  Despite  the  constraints,  the  controller  is  still  able  to  guide  the 
spacecraft  to  the  desired  attitude. 
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(a)  (b) 

Figure  5,  MFC  Response  of  Spacecraft  Controlled  by  two  Internal  Moments  Produced  by 
RWs  to  (a)  Euler  Angles,  (e)  Control  Accelerations 


5.0  CONCLUSION 

This  work  demonstrated  several  approaehes  to  the  attitude  eontrol  of  the  underaetuated  spaeeeraft. 
The  first  method  utilizes  SRP  to  regain  linear  eontrollability  of  an  underaetuated  spaeeeraft.  With 
linear  controllability  restored,  conventional  control  techniques  such  as  LQ  control  can  be  applied. 
The  second  method  demonstrates  the  ability  to  construct  sets  of  initial  conditions  for  which  if  a 
spacecraft  starts  within,  it  will  satisfy  pointing  and  state  constraints  regardless  of  the  disturbance 
torque  exerted  on  it.  The  third  method  demonstrates  the  use  of  hybrid  and  MPC  schemes  for  the 
underaetuated  spacecraft  attitude  control.  Both  approaches  deal  with  nonzero  angular 
momentum,  which  is  the  case  usually  not  considered  in  the  current  literature,  yet  has 
important  consequences  with  regards  to  the  operational  life  of  the  spacecraft.  Both  controllers 
are  able  to  successfully  guide  to  the  attitude  equilibrium. 

Future  work  will  be  done  on  the  development  of  the  drift-sets  approach.  Constructing  these  sets 
can  be  computationally  intensive,  so  the  use  of  other  methods  such  as  those  based  on  zonotopes  is 
of  great  interest.  In  addition,  utilizing  these  drift  sets  as  a  means  to  construct  a  control  law  is  the 
next  step  of  the  evolution  of  the  work.  The  MPC  approach  for  the  zero  angular  momentum  case 
will  also  be  studied  more  in  depth,  as  a  nonlinear  controller  that  provides  for  rapid  implementation 
will  need  to  be  developed. 
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LIST  OF  SYMBOLS,  ABBREVIATIONS,  AND  ACRONYMS 


FUSE:  Far  Ultraviolet  Spectroscopic  Explorer 

EQ.  Einear  Quadratic 

RW:  Reaction  Wheels 

SRP:  Solar  Radiation  Pressure 

MPC:  Model  Predictive  Cntroller 

0:  Roll  Euler  angle 

9:  Pitch  Euler  angle 

0:  Yaw  Euler  angle 

0:  Vector  of  stacked  Euler  angles 

M(0):  Euler  angles  kinematics  matrix 

o):  Angular  velocity 

J :  Inertia  matrix  of  the  spacecraft  bus 

J„:  Inertia  of  the  RW  about  the  spin  axis 

W :  Orientation  of  the  RW  spin  axes  with  respect  to  the  spacecraft  bus 
v:  RW  speeds 

Mgxf-:  External  moments  acting  on  the  spacecraft 

OgQ :  Direction  cosine  matrix  from  the  inertial  frame  to  the  spacecraft  frame 
H:  Inertial  angular  momentum  vector 
Pj :  SRP  pressure  exerted  on  panel  j 
a:  SRP  constant 
P:  SRP  constant 

Ujij:  Normal  to  panel  j  point  outward 
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Unit  vector  from  the  center-of-mass  to  the  sun 
Tj :  SRP  torque  exerted  by  panel  j 

Tjo :  Vector  from  center-of-mass  to  the  center  of  panel  of  j 

Tco-  Vector  from  the  center-of-pressure  to  the  center-of-mass 

Aj :  Area  of  panel  j 

x(A:):  State  at  time  instant  k 

u{ky.  Control  at  time  instant /c 

w{ky.  Bounded  disturbance  at  time  instant  k 

A:  Discrete  state  matrix 

B  :  Discrete  control  matrix 

U :  Set  of  acceptable  control 

W :  Bound  on  disturbance 

Hi,  hi'.  Linear  equability  definitions 

n:  Number  of  linear  equalities 

0]^:  Set  of  all  initial  conditions  for  which,  if  the  spacecraft  starts  with,  satisfies  constraints  at  time 
k 

^sun-  Solar  Flux 

Cdiff'-  Diffusion  coefficient 

c  ;  Speed  of  light 
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